Slave-boson approach to the metallic stripe phases with large unit cells 
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Using a rotationally invariant version of the slave-boson approach in spin space we analyze the 
stability of stripe phases with large unit cells in the two-dimensional Hubbard model. This approach 
allows one to treat strong electron correlations in the stripe phases relevant in the low doping regime, 
and gives results representative of the thermodynamic limit. Thereby we resolve the longstanding 
controversy concerning the role played by the kinetic energy in stripe phases. While the transverse 
hopping across the domain walls yields the largest kinetic energy gain in the case of the insulating 
stripes with one hole per site, the holes propagating along the domain walls stabilize the metallic 
vertical stripes with one hole per two sites, as observed in the cuprates. We also show that a finite 
next-nearest neighbor hopping t' can tip the energy balance between the filled diagonal and half- 
filled vertical stripes, which might explain a change in the spatial orientation of stripes observed in 
the high T c cuprates at the doping x ~ 1/16. 

PACS numbers: 74.72.-h, 71.10.Fd, 71.45.Lr, 75.10.Lp 



INTRODUCTION 



The experimentally established existence of nanoscale 
charge order phenomena in transition metal oxides con- 
tinues to attract much attention. In particular, phase 
separation, manifesting itself in formation of nonmag- 
netic one-dimensional (ID) domain walls which sepa- 
rate the antiferromagnetic (AF) domains with opposite 
phases, has been proposed to be responsible for unusual 
properties of the high T c copper oxide superconductors. 1 
Moreover, taking into account that magnetic excitation 
spectra have been found to be similar among various high 
T c cuprates, the knowledge of the spatial distribution of 
holes in the CuC"2 planes might be a first step towards 
microscopic understanding mechanism of the supercon- 
ductivity itselfp^ Indeed, the most prominent feature of 
the spin excitations in YBa2Cu306+a (YBCO) is the res- 
onance that occurs at the AF wavevector Q = (n, ir) and 
at the energy of 41 meV in the case of optimal doping 
while for lower energies incommensurate (IC) magnetic 
peaks emerged In contrast, for a long time it was thought 
that La2- a: Sr a ;Cu04 (LSCO) does not show the resonance 
and that the corresponding IC fluctuations are disper- 
sionless. However, recent high-resolution neutron scat- 
tering studies on optimally doped sample have revealed 
that in fact the magnetic excitations are dispersive and 
quite similar to those of YBCOi^ Moreover, the over- 
all evolution of the magnetic scattering with energy in 



La2- a; Ba 2 ;Cu04 (LBCO) also resembles the one found in 
YBCO£ Altogether, the common low-energy excitations 
imply that the spin dynamics in doped copper oxide su- 
perconductors has the same origin. Remarkably, the ob- 
served spectra might be consistently explained in terms of 
fluctuating stripes suggesting that the stripe instability is 
a universal phenomenon of all cuprates and thus it might 
be indeed regarded as relevant for the superconductivity^ 
Historically, the first compelling evidence for both 
magnetic and charge order in the cuprates, was 
accomplished in a neodymium codoped compound 
Lai.6-xNdo.4Sra;Cu04 (Nd-LSCO). Indeed, around hole 
doping x — 1/8, Tranquada et alJ found that elastic 
neutron scattering is characterized by magnetic peaks at 
the wave vectors Q = 7r(l ± 2e, 1) and Q = tt(1, 1 ± 2e). 
Such positions of the peaks corresponds to equally prob- 
able modulations along one of two equivalent lattice di- 
rections x and y, to which we refer as horizontal and 
vertical stripes, respectively. Moreover, inspired by the 
pioneering Hartree-Fock (HF) studies suggesting that the 
staggered magnetization undergoes a phase shift of n at 
the charge domain wall (DW)r the authors found ad- 
ditional Bragg peaks displaced symmetrically by ±4ire 
around the T = (0, 0) point, precisely at the position 
expected for charge modulation. Recently, the same ge- 
ometry for static IC magnetic and charge peaks has been 
reported^, in LBCO with x = 1/8. In fact, neutrons do 
not detect charge directly but instead they are sensitive 
to nuclear displacements induced by the charge modula- 
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tion. Therefore, a very recent resonant soft X-ray scat- 
tering study of the charge order in LB CO at the same 
doping is of a particular great value being a more di- 
rect evidence of charge modulation. These studies have 
yielded an estimate for the number of holes per DW to be 
0.59, very close to the expected from neutron scattering 
experiments value 0.5. It corresponds to the so-called 
half-filled stripes, as they are characterized by the fill- 
ing of one doped hole per two atoms along the DWs. 10 
Remarkably, at this particular doping, both compounds 
exhibit a deep minimum in the doping dependence of T c , 
suggesting a strong competition between the static stripe 
order and the superconductivity^ 

Conversely, the large drop in T c at x — 1/8 is not ob- 
served in LSCO, even though inelastic neutron scattering 
experiments in the superconducting regime of x > 0.06 
have revealed the presence of magnetic peaks at the same 
IC positions as in Nd-LSCO and LBCO. Moreover, Ya- 
mada et al. 1 - established a remarkably simple relation 
e ~ x for 0.06 < x < 0.125 with a lock- in effect at 
e ~ 1/8 for larger x. This means that increasing dop- 
ing reduces the distance d = l/2x between the half-filled 
stripes. Hence, it appears that such a stripe order is com- 
patible with the superconductivity, which is particularly 
pronounced when the spin correlations remain purely dy- 
namic. However, the superconducting state may also co- 
exist with static stripe order^ but then the value of 
T c is strongly reduced. This conjecture is strongly sup- 
ported by inelastic neutron scattering experiments on 
YBCO, which have also established the presence of IC 
vertical/horizontal spin fluctuations throughout the en- 
tire superconducting regime of YBCOi 1 ^ Moreover, for 
S — 0.35, apart from spin fluctuations, Mook et al^i 
have found an IC charge order consistent with the verti- 
cal/horizontal stripes. 

In contrast, in the insulating spin-glass regime of LSCO 
x < 0.06, quasielastic neutron scattering experiments 
with the main weight at zero frequency demonstrate 
that IC magnetic peaks are located at the wave vector 
Q = 7r(l ± V2e, 1 ± \/2e)iA£ This phenomenon has of- 
ten been interpreted as the existence of static diagonal 
stripes, even though no signatures of any charge mod- 
ulation were observed. In spite of the change in spin 
modulation from the diagonal to vertical/horizontal one 
at the doping x ~ 0.06, e follows x reasonably well over 
the entire low doping range 0.03 < x < 0.125. More- 
over, the same type of the diagonal IC spin modulation 
has been discovered in Nd-LSCO at x — 0.05, 17 indicat- 
ing that while the vertical/horizontal spin-density mod- 
ulation, accompanied by the charge modulation, might 
be regarded as a generic property of the superconduct- 
ing regime, the diagonal orientation of the magnetic 
domains is common for the lightly Sr-doped insulating 
systems. Finally, in the narrow range of low doping 
0.02 < x < 0.024, IC magnetic peaks in LSCO are ob- 
served at the wave vector Q = 7r(l ± e, 1 ± e)^ just as 
in the insulating La2- a; Sr 2 ;Ni04 (LSNO) compound^ In 
this case, the spacing between stripes is equal to d = 1/x; 



such structures with the filling of one doped hole per one 
atom in the DW correspond to the so-called filled stripes. 

The intriguing competition between vertical and diag- 
onal stripes has been already noticed in the early mean- 
field studies^ Indeed, it has been found that the vertical 
stripes were favored for U/t < 4 whereas the diagonal 
DWs were formed at a stronger Coulomb repulsion. Un- 
fortunately, these studies predicted the filled stripes in 
the ground state. A general feature of such instability 
is a gap/pseudogap precisely on the Fermi surface. As a 
result, charge transport is not possible in idealized filled 
stripes. In contrast, the ground state energy half-filled 
stripes observed in the cuprates, have been found in a 
few methods which go beyond the Hartree-Fock approx- 
imation (HFA), such as: density matrix renormalization 
group^fi variational local ansatz approximation j2I Exact 
Diagonalization (ED) of finite clusters^ analytical ap- 
proaches based on variational trial wave function within 
the string picture ; 23 ' 24 dynamical mean field theory^ 
cluster perturbation theory^ and Quantum Monte Carlo 
(QMC) simulations^ indicating the crucial role of local 
electron correlations in stabilizing these structures. In- 
deed, while the HF method is well suited to compare the 
energies of different types of magnetically ordered phases 
in the limit of large on-site Coulomb interaction U — > oo, 
where it gives the same asymptotic behavior as the SBA 
or the Gutzwiller ansatz ^ it is crucial to use a better 
approach than the HFA at intermediate U, particularly 
when magnetic polarization of some sites is missing. 

A good starting point for a proper approximate treat- 
ment of strong correlations is the representation of 
the Hubbard model in terms of auxiliary fermions and 
bosons? 2 ^ The slave-boson approximation (SBA) has 
been applied successfully to a whole range of problems 
and is known to provide a realistic mean-field descrip- 
tion of strongly correlated systems. It is quite encourag- 
ing that the ground state phase diagram obtained using 
spin-rotation-invariant (SRI) slave-boson representation 
of the two-dimensional (2D) Hubbard model with homo- 
geneous spiral, AF order, ferromagnetic (FM), and para- 
magnetic (PM) phases shows a good agreement both with 
QMC simulations and the ED method^ The slave-boson 
(SB) method was also used to investigate magnetic and 
charge correlations of the t-t'-XJ model, 31 the ground state 
of the Anderson lattice modelj^ and systems with orbital 
degeneracy^ Moreover, the unrestricted slave-boson for- 
malism has turned out to be a powerful tool in description 
of inhomogeneous states, even though in the absence of a 
finite long-range Coulomb repulsion, completely filled di- 
agonal stripes (FDS) were found to be more stable than 
the half-filled vertical stripes (HVS), suggesting that a 
pure Hubbard model might be insufficient to capture the 
physics in the cuprates i 34 ' 3 ^ 

Another possible extension could involve an inclusion 
of the next-nearest neighbor hopping term t' . There are 
several experimental and theoretical studies suggesting 
the presence of a finite t' in the cuprates. Indeed, topol- 
ogy of the Fermi surface seen by angle-resolved photoe- 
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mission spectroscopy 3 - 6 and the asymmetry of the phase 
diagrams of the hole- and electron-doped cuprates can be 
understood only by introducing f'i2X It also offers an ex- 
planation for the variation of T c among different families 
of hole-doped cuprates j 38 ' 39 Moreover, ED studies have 
shown that while the d-wave superconductivity correla- 
tion is slightly suppressed by f in underdoped regions, 
it is substantially enhanced in the optimally doped and 
overdoped regions, indicating that t' is of great impor- 
tance for the pairing instability^ 

Finally, Seibold and Lorenzana have found^i that a fi- 
nite t' strongly affects the optimal filling of the vertical 
stripes and consequently favors partially filled configura- 
tions in the SBA. Unfortunately, restrictions due to the 
cluster size, have not let the authors to reach an unam- 
biguous conclusion concerning a crossover towards the 
FDS at low doping. Motivated by this result, we have 
investigated recently the relative stability between the 
latter and the HVS with the symmetry lowered by a pe- 
riod quadrupling due to on- wall spin-density wave. Using 
the self-consistent HFA, one finds that a negative ratio 
of the next-to nearest neighbor hopping (f'/f < 0) expels 
holes from the AF domains and reinforces the stripe or- 
der. Therefore, the half-filled stripes not only accommo- 
date holes but also redistribute them so that the kinetic 
energy is gained, and these stripes take over in the regime 
of t'/t ~ -0.3 appropriate for YBC042 

The purpose of this paper it to investigate the stability 
of stripe phases at low doping and to resolve the long- 
standing controversy concerning the role of the kinetic 
energy in the formation of stripes. Domain walls can 
be viewed as topological defects in the AF phase, which 
confine the holes to ID regions, while a crossover to a 2D 
system occurs when the stripes melt4£ This leads in a 
natural way to the picture of ID composite excitations^ 
which contribute to the kinetic energy by the hole motion 
along DWs. In contrast, in the early HF studies it was 
demonstrated that the stripes are stabilized by the ki- 
netic energy increments due to the bonds connecting the 
domain wall with the neighboring sites of AF domains. 44 

The rest of the paper is as follows. In Sec.[H]we shortly 
review the SRI slave-boson representation of the Hub- 
bard model and introduce reciprocal space formalism, 
based on the stripe periodicity, which provides a possibil- 
ity to perform calculations on larger (128x128) clusters 
than those studied in Ref. \^ Therefore, our approach 
eliminates to a large extent the role of finite-size effects, 
as we show it in Sec. Ill PI Further, in Sec. lIIII we investi- 
gate the stability of the already mentioned FDS and HVS 
as well as of the filled vertical stripes (FVS) and half-filled 
diagonal stripes (HDS) at the two representative doping 
levels x = 1/8 and x — 1/16. For stable structures we 
compare the charge, spin, and double occupancy profiles 
obtained in the SBA with the ones obtained in the HFA. 
We also analyze in detail kinetic and magnetic energy 
contribution as well as discuss the shape of the density 
of states (DOS). The effect of a finite next-nearest neigh- 
bor hopping t' on the stripe stability and induced by it 



conspicuous changes in the band structure are presented 
in Sec. IIVI Finally, in Sec. 0we present a short sum- 
mary and general conclusions. 

II. FORMALISM 
A. Rotationally invariant slave-boson method 

In this paper we study the 2D Hubbard model, 

H = -Y^Ujc\ a c ja + U^n i ^n ii , (1) 

ij,<r i 

where the electron hopping t t j is t on the bonds con- 
necting the nearest neighbor sites (ij) and t' between 
next-nearest neighbor sites, while U stands for the 
on-site Coulomb interaction. Our study is based on 
the spin-rotation invariant formulation of the Kotliar 
and Ruckenstein^S. SB representation of the Hubbard 
model) 4 *** 4 ^ in which one enlarges the Hilbert space of 
the Hamiltonian Q by introducing auxiliary boson op- 
erators: and rfj, as well as a scalar p iQ and a vector 
Pi = (Pii)Pi2iPi3) bosonic field at each site. The for- 
mer operators act as projection operators on empty and 
doubly occupied sites, while the latter describe the spin 
and charge degrees of freedom in the singly occupied sub- 
space. 

The main advantage of such a representation is that the 
actual electron configuration is controlled by the bosons 
and one can thus write the Hubbard interaction as a 
bosonic occupation number operator. In contrast, the 
operator for the kinetic energy becomes more involved 
since the motion of a physical electron changes boson oc- 
cupation numbers on both lattice sites involved in the 
hopping process. Hence, in terms of the above SB oper- 
ators, the Hubbard model Q takes the form, 

ij a a' a" i 

(2) 

where are 2x2 matrices in spin space which depend on 
the actual configuration of the boson fields, as explained 
in Ref. 53 These matrices are constructed in such a way 
that the correct mean-field result in the noninteracting 
limit (at U = 0) is recovered. Additionally, the SB oper- 
ators have to fulfill the following constraints at each site, 

24^+E4^=E/*U- ( 4 ) 

PoiPi + p!po* ~ «P! X P* = E T ™' fi°' fi°' ( 5 ) 



4 



They are enforced by corresponding Lagrange multipliers 
in the action, and one obtains a path integral representa- 
tion of the Hubbard model. Owing to the gauge symme- 
try group of the action, the phase of five of the six bosonic 
fields can be gauged away4& These fields can be handled 
in the Cartesian gauge, or in the radial gauge4&In the 
following we handle the action in the Cartesian gauge at 
saddle-point level (in the SB A). This procedure is exact 
in the large degeneracy limit Comparison of ground 
state energies with data resulting from numerical simu- 
lations show very good agreement £M£L 

Since homogeneous AF phases turned out to be un- 
stable with respect to spiral phasespS, one may antici- 
pate that a similar effect will take place in stripe phases. 
Therefore, the present rotationally invariant formalism 
opens a possibility of studying also more general struc- 
tures which we postpone to future studies. Here we will 
limit ourselves to noncanted states which are locally sta- 
ble, as checked in a few representative cases. 



B. Reciprocal space approach for stripe phases 

In order to obtain unbiased results for stripe phases one 
should carry out calculations on large clusters. So far, 
such calculations were usually performed in real space 
to address directly inequivalent sites in the stripe unit 
cell. However, an efficient approach may be constructed 
by making use of the periodicity of a stripe phase which 
allows one to cover the entire lattice by small unit cells. 
Indeed, each position vector Ri of an arbitrary atom can 
be decomposed into a sum involving certain periodicity 
vectors g ± and g 2 , and a vector 8m labeling all inequiv- 
alent sites i within the unit cell: 



(6) 



with integer {n, m}. Accordingly the covering of the Bril- 
louin zone spanned by the wave vector k is achieved by 
the decomposition: 



k = K + G, 



(7) 



with G belonging to a set of wave vectors representative 
of the periodicity of the structure of interest £SL For a 
stripe phase with vertical domain walls separated by d = 
4 lattice spacings, as shown in Fig. QJa), the magnetic 
unit cell consists of a row of eight atoms perpendicular 
to the stripes. One immediately finds two periodicity 
vectors g x — (4, l) and g 2 = (0, 2J . In this case the set 
of vectors G satisfying expjig^ 2 ) • G} = 1 is given by 



G« = |(2/,0), 



G 



(') 



?(« + M). 



(8) 



with < I < 3. 

A second example is a diagonal stripe phase depicted 
in Fig. Gib). Again, considering the periodicity of this 
phase, one finds the magnetic unit cell consisting of a 
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FIG. 1: Symmetry properties of stripe phases: (a) the vertical 
stripe phase — its unit cell with two periodicity vectors: g 1 = 
(4, 1) and g 2 = (0, 2) ; (b) the diagonal stripe phase — its 
unit cell with two periodicity vectors: g 1 = (1, 1) and g 2 = 
(4,-4). 



row of eight atoms along the x-direction. However, the 
periodicity vectors are now given by g x — (1,1) and g 2 = 
(4, —4). Accordingly one finds for G: 

2tt 



Gi = (0,0), 



i(0 



-(1,8-0, 



(9) 



with 1 < I < 7. 

In both cases the vectors K in Eq. Q are spanning 
the rectangle [0, 2ir/8] x [0, 2tt]. The extension to other 
stripe phases is straightforward. Using the above pro- 
cedure, the fermionic contribution to the action in the 
momentum representation is given by: 



S f = 



with, 



f dT E E/K+G,.M^G'(K,r)/ K+G> ,(r), 

J ° KGG' era' 

(10) 



M^ G ,(K,r) = {d T -n)5 a „,5 GGI + 1 ^ e -i(G-G')8 M 

• aa' 

:ts M ,s M+ se« K+G '^ (z SM+s (r)4jr)) T 



s 



(11) 
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Here L is the number of atoms in the unit cell so that the 
ratio N/L gives the number of unit cells needed to cover 
a whole cluster with N sites, d T is an imaginary time 
derivative, and [i is the chemical potential. The quanti- 
ties /3o and (3 are the Lagrange multipliers enforcing the 
constraints given in Eq. £Q and Eq. ©, respectively. 

Hence, we have reduced the fermionic matrix M££, 
down to decoupled blocks labeled by K. This procedure 
gives a considerable time gain during numerical calcu- 
lations. While the usual workload in a diagonalization 
algorithm is ~ TV 3 (see Ref. 02), the number of opera- 
tions needed for diagonalization of N/ L smaller matrices 
is only ~ =r • L 3 = NL 2 . This means that the sym- 
metry reduction makes a tremendous simplification and 
the time needed for numerical calculations is reduced by 
a factor (N/L) 2 as compared to a straightforward 'brute 
force' diagonalization of an TV x N matrix, which is neces- 
sary when the computations are performed in real space. 



C. Mean-field approximation 

In Sees. Ill and IV we present solutions for the stripe 
phases obtained within the SB mean-field approximation 
in which one replaces the Bose fields and the Lagrange 
multipliers by their time-independent averages. Hence, 
the free energy found at the temperature T follows as^ 30 

I 

- ftoMi +Pf + 2d?) + Udf - 2/3, • PlPol ] 



-/3e q 



ft 



(12) 



with (3 = 1/fcsT, e qCT standing for the eigenmodes of 
the fermionic matrix 11 H , and N e i corresponding to the 
total number of particles. The equilibrium values of the 
classical field amplitudes are determined from the saddle- 
point equations in which the partial derivatives are taken 
with respect to three Lagrange multipliers enforcing the 
constraints 0~0> as well as to the four SB fields poi, pi, 
e/, and di, for each inequivalent site I = (lx,ly) within the 
elementary stripe unit cell. The saddle-point equations 
are solved using the Powell hybrid method. 



D. 



Finite-size effects 



In order to estimate to what extent the calculation 
method eliminates finite-size effects, we examine the size 
dependence of the free energy at low temperature ftt = 
1000 by considering stripe phases on clusters of increasing 
size: 8 x 8, 16 x 16, 24 x 24, 32 x 32, 64 x 64, and 128 x 128. 
All systems are described by the Hubbard model with 
U = 12t and x = 1/8. The finite-size scaling of the free 
energy obtained for the HVS obtained within the SBA is 
shown in Fig. EI a). As it could be expected, finite size 
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FIG. 2: Finite-size scaling of free energy F for: (a) HVS, and 
(b) HDS, as obtained within the SBA for the Hubbard model 
at doping x = 1/8. Parameters: U = 12t and fit = 1000. 



effects are particularly severe for small 8x8 and 16 x 16 
clusters but a further increase of the system size results 
in a gradual saturation of the free energy. The systematic 
errors due to these effects are diminished to the order of 
10~ 5 i when the calculations are performed on clusters 
with 10 4 atoms. 



The above behavior is universal, and a similar tendency 
appears also in a system with the HDS [see Fig. E^b)]. 
Note, however, that in the case of the vertical stripes 
increasing cluster size lowers the total energy. In con- 
trast, the energy increases when the number of atoms in 
a cluster with the diagonal stripes gets larger, suggesting 
that small cluster calculations underestimate (overesti- 
mate) the free energy of the vertical (diagonal) stripe 
phase, respectively, and may thus lead to systematic er- 
rors and inconclusive results in some cases. Even though 
the free energy F saturates quickly with an increasing 
number of sites in both cases, one should keep the above 
feature in mind, especially when it comes to investigate 
the relative stability of vertical and diagonal phases as 
the inclusion of quantum fluctuations could substantially 
strengthen this effect. Here we would like to point out 
another virtue of the present studies performed at tem- 
perature ftt = 1000. When the temperature is so low, the 
entropy contribution to the free energy is almost entirely 
suppressed as F ~ (H), so that one may analyze only 
the internal energy for different phases, as done below. 
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III. ELECTRON CORRELATION EFFECTS IN 
STRIPE PHASES 

A. Physical characterization of stripe phases 

Let us now examine the role of strong electron corre- 
lations in stabilizing stripe phases of various type. Even 
though the HFA and SBA yield quantitatively very sim- 
ilar results at half-filling, as recently reviewed by Kor- 
bel et al.jH we show below that they differ markedly in 
the stripe phases. Indeed, since the HF wave function 
is a Slater determinant of single particle states, it does 
not provide enough variational freedom to implement the 
many-body processes that are relevant for the behavior of 
interacting fermions. Therefore, one should expect that 
the inclusion of the correlation effects modifies consid- 
erably the distribution of charge and spin density in a 
stripe phase, especially around nonmagnetic DWs where 
the correlation corrections are largei^I 

To illustrate this point, we investigate the charge and 
magnetization distribution in the SBA and compare them 
with the reference states obtained in the HFA. To de- 
scribe the charge distribution, we introduce the local hole 
density for each inequivalent site I = (l x ,l y ) in the stripe 
unit cell, 

n h (l x ) = 1 - J2(fiM> ( 13 ) 

(7 

where the average electron density, 

(flU = E *l(q)^,(qM^), (14) 

q 

is given by the real space eigenvectors ^^(q) of the ef- 
fective Hamiltonian llOII . and n(£ qcr ) is the Fermi-Dirac 
distribution. Further, magnetic domain structures are 
best described by the modulated magnetization density, 

sm = E<//^-^'>> ( 15 ) 

with a site dependent phase factor (— 1)'*+'* compensat- 
ing the modulation of the staggered magnetization den- 
sity within a single AF domain, as well as by double 
occupancy, 

D(l x ) = (m T n u >. (16) 

B. Stripes at low doping x = 1/16 

We begin with the detailed analysis of the charge and 
magnetization distribution in the stripe phases at low 
doping x = 1/16, and compare the results of the SB treat- 
ment with those obtained in the HFA. The corresponding 
SB profiles of the FVS (left) and FDS (right) found in the 
Hubbard model with U = 6t at the doping x = 1/16 are 
shown in Fig. (filled circles). In this case nonmagnetic 



DWs are separated by d = 16 lattice constants so that 
the charge (magnetic) unit cell contains sixteen (thirty 
two) atoms, res_pectively. In agreement with the calcu- 
lations of Ref. 1341 we note that the hole density n% at 
nonmagnetic DWs is reduced nearly twice in the SBA 
as compared to the corresponding HFA value (open cir- 
cles), regardless of the stripe direction. Moreover, the HF 
stripe phases are characterized by the enhanced spin po- 
larization of the AF domains. Such a strong modification 
follows directly from the fact that in the HFA there are 
only two straightforward possibilities which allow one to 
diminish the energy of the stripe phase due to the on-site 
Coulomb repulsion, 

Eu = U~J2D(l x ). (17) 

The simplest way of keeping apart electrons with the op- 
posite spins is by creating strong spin polarization. As 
a consequence, a well known feature of the HFA is that 
it overestimates by far the tendency towards symmetry 
broken states. Another way of reducing Exj which may 
be realized for inhomogeneous charge distribution, for 
instance in stripe phases, is to suppress locally the to- 
tal electron density at the sites with small or vanishing 
magnetization. 

In contrast, the SB approach implements local correla- 
tions by offering an important mechanism to optimize the 
on-site interaction by an additional variational parame- 
ter, i.e., the boson field cf"j at each site. In fact, the largest 
correction of the HF value is obtained at the DW unpo- 
larized sites, where the double occupancy D, Eq. ltl6)l . 
shows distinct minima in both phases (see Fig. which 
allows one to optimize Ey even without a great reduction 
of the actual electron density. Simultaneously, a large 
value of D within the AF domains suppresses partially 
the spin polarization of the atoms separating DWs, which 
enables more intersite excitations and leads to a more fa- 
vorable kinetic energy gain. Taken together, these two 
effects are jointly responsible for smoother charge and 
spin density profiles with respect to the ones found in 
the HFA. Notice that both approximations yield narrow 
diagonal stripes, revealing their more localized character 
as compared to the vertical ones and hence a more favor- 
able average on-site energy, which stabilizes them in the 
strong coupling regime, as discussed below. 

It is quite remarkable that the present calculations per- 
formed on large clusters yield also locally stable half-filled 
stripes both in the HFA and in the SBA without any ne- 
cessity of quadrupling of the period along the stripes by 
an additional on-wall spin-density wave4^ At x = 1/16, 
the condition of having one doped hole per two DW 
atoms requires the charge (magnetic) unit cell containing 
eight (sixteen) atoms, respectively. The obtained profiles 
of the HVS (left) and HDS (right) are shown in Fig. H 
Again, one finds a smaller SB charge modulation and a 
stronger spin polarization of the AF domains in the HFA. 
Note, however, that contrary to the case with the filled 
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FIG. 3: (Color online) Hole and magnetization distribution 
in filled stripe phases at low doping x — 1/16: local hole 
n h (Z:r) density (top), magnetization S^ilx) density (second 
row), and double occupancy D(l x ) (bottom), for the phase 
with the FVS (left) and FDS (right), found in the Hubbard 
model on a 128 x 128 cluster at (3t = 1000. Open (filled) 
circles show the HFA (SBA) results for U = 6t, respectively, 
while filled squares denotes the SBA data for U = 12t. 



DWs, a narrower charge and spin profile of the vertical 
stripe with respect to the diagonal one is apparent. 

At strong coupling with U = 12t (filled squares), one 
finds for the most stable phases that the double occu- 
pancy profile is interpolating between D ~ At 2 /U 2 deep 
in the magnetic domains (the value expected in an AF 
phase at half-fillingS) , and the value expected in the PM 
phase at the actual hole density realized at the domain 
wall sites, as shown in Table[I] Thus the present approach 
is flexible enough to implement the compromise between 
optimized superexchange in the magnetic domains, and 
optimized (with respect to the local density) double oc- 
cupancy in the domain walls, with, however, the notable 
exception of the HDS. In that case this optimization can- 
not be achieved, and instead one gets a more spread hole 
distribution. The same tendency is also realized at inter- 
mediate coupling. 



FIG. 4: (Color online) The same as in Fig.|2]but for the half- 
filled stripe phases: HVS (left) and HDS (right). Open (filled) 
circles show the HF (SB) results for U = 6t, respectively, while 
filled squares denotes the SB data for U = 12t. High values 
of UD/t > 0.6 obtained for the DW atoms in the HFA are 
indicated by arrows. 



Experimentally, two stripe phases have been found: (i) 
filled diagonal stripe phase, and (ii) half-filled vertical 
one. Both of them have been observed in LSCO; the 
former was found around x = 0.02p£ whereas the latter 



TABLE I: Hole density iih(DW) and double occupancy 
D(DW) at domain-wall sites, and the corresponding dou- 
ble occupancy in the PM phase Dpm, as found for various 
stripe phases in the Hubbard model on an 128 x 128 cluster 
at x = 1/16 doping within the SBA. Parameters: U = 12t, 
[3t = 1000. 



phase 


%(DW) 


JJ(DW) 


-Dpm 


HDS 


0.153 


0.016 


0.019 


HVS 


0.221 


0.018 


0.018 


FVS 


0.256 


0.017 


0.017 


FDS 


0.267 


0.017 


0.016 
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— at a higher 0.06 < x < 0.125 doping levels In the 
present calculations these phases were reproduced as two 
possible ground states of the doped Hubbard model. The 
crossover found experimentally indicates that the filling 
and orientation of DWs are indeed closely related in the 
cuprates. In other words, half- filled stripes tend to be 
aligned vertically /horizontally, whereas filled ones gain 
more energy by the diagonal arrangement. 

The above important property can also be deduced 
from the free energies F calculated in the SBA from 
Eq. ltl2ll for all four stripe phases, given in Table [H] 
For completeness, the energies of stripe phases are com- 
pared with those of the uniform AF phase. Note that 
the sequence of stripe phases, ordered according to their 
decreasing energy, is precisely the same in both methods 
and one does find that the half-filled vertical stripe phase 
is preferred over its diagonal counterpart, and also that 
the filled diagonal stripe structure is promoted over the 
corresponding vertical stripe phase. However, noticeable 
discrepancies emerge when it comes to compare the on- 
site Coulomb energy Ejj and the nearest neighbor kinetic 
energy determined in SBA, 

^=-4eem)L<4/^>, as) 

{1,1') aa' 

with the ones obtained in the HFA. Most importantly, 
the inclusion of electron correlations reduces Ejj of the 
PM phase by a factor two without a strong suppression 
of E t , which improves significantly its free energy. There- 
fore, as the average hole density at the half-filled DWs is 
noticeably smaller than at the filled ones (c/. Figs. 
andQJ, it should also help to stabilize half-filled stripes, 
because one gains more correlation energy when the non- 
magnetic atoms are close to half-filling. These features 
are clearly seen in Table [H] By comparing the repulsion 
energy, one finds that Eu of both half-filled stripe phases 
is reduced in the SBA, while Eu of the filled ones is even 
slightly enhanced. As a result, for moderate coupling, 
the structure with the HVS in the SBA is characterized 
by the smallest net double occupancy and hence by the 
most favorable Eu- Nevertheless, in the SB ground state 
one recovers the phase with the FDS, as a small energy 
loss due to the increased on-site energy is easily over- 
compensated by the kinetic energy gain. The FDS are 
recovered at strong coupling as well, where the HVS are 
not even favored by Eu. 

Finally, it should be emphasized that, in contrast to 
the HFA, the most favorable SB phase in the moderate 
coupling regime U = 6t for the intersite kinetic energy 
E t is the AF structure (see Table ITTf) . However, owing 
to a strong suppression of the double occupancy at non- 
magnetic DWs, the existence of stripes opens a possibil- 
ity for a charge redistribution which optimizes Eu even 
better. Consequently, in the strong coupling regime with 
U = 12i all the stripe phases are favored over the uniform 
AF order. 

In order to identify in detail processes leading to the 
stripe stabilization, we determined the expectation values 



of the bond hopping terms, E'f(l x ) and E^(l x ), along the 
x and y direction, respectively, 

E*{l x ) = -*E^(M</ ( b H ),./fc +1 ^),^ (19) 

era' 

where, 

C'fe) = (^ +1 ,j H )4„w)L' (21) 

c<(y = (* ( w* + o4^))l' (22) 

are the corresponding band narrowing factors. The vari- 
ation of the kinetic energy gain across the unit cell of 
the phase with the FVS is shown in Figs. EI a) and[^b). 
Deep in the magnetic domains, for large U, the kinetic 
energy in both directions is given by —At 2 /U, as expected 
from the large U expansion^ For the DW sites we obtain 
that Ef closely reproduces one half of the kinetic energy 
in the PM phase at the actual hole doping found at these 
sites. In the two cases shown in Fig. the kinetic en- 
ergy is weakly influenced by the involved band structure 
describing the stripe phases, and its profile mostly inter- 
polates between the two expected limits. For the FVS its 
most significant effect is seen in Ef close to the DWs, 
where most of the kinetic energy gain is realized. The 
same trends are also found for smaller U. 

As expected, profile of Ef(l x ) and E^(l, x ) depends on 
whether we include local electron correlations or not, as 



TABLE II: Comparison of the site-normalized free energy 
F, on-site Coulomb energy Eu, and kinetic energy E t , as 
obtained for various phases in the Hubbard model within the 
HFA (for U = 6t) and SBA (for U = 6t and 12t) on an 128 x 
128 cluster at doping x — 1/16 and (5t = 1000. For clarity, 
the phases are listed according to their decreasing energy from 
top to bottom. 



method 


U/t 


phase 


Eu/t 


Et/t 


F/t 


HFA 


6 


HDS 


0.5185 


-1.1080 


-0.5895 






AF 


0.4928 


-1.1198 


-0.6270 






HVS 


0.4892 


-1.1242 


-0.6350 






FVS 


0.4690 


-1.1218 


-0.6528 






FDS 


0.4586 


-1.1154 


-0.6568 


SBA 


6 


HDS 


0.4796 


-1.1708 


-0.6912 






AF 


0.4943 


-1.2009 


-0.7066 






HVS 


0.4663 


-1.1809 


-0.7146 






FVS 


0.4709 


-1.1929 


-0.7220 






FDS 


0.4671 


-1.1902 


-0.7231 


SBA 


12 


AF 


0.3181 


-0.7355 


-0.4174 






HDS 


0.2673 


-0.7023 


-0.4350 






HVS 


0.2860 


-0.7356 


-0.4496 






FVS 


0.2930 


-0.7455 


-0.4525 






FDS 


0.2903 


-0.7464 


-0.4561 
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FIG. 5: (Color online) Kinetic energy projected on the bonds 
in the x (y) direction, Ef(l x ) and Ef(l x ), as found in the phase 
with either the FVS (a,b) or HVS (c,d), shown in Figs. H 
and El Open (filled) circles show the HF (SB) results for 
U = 6t, respectively, while filled squares denote the SB data 
for U = 12t. Panels (e) and (f) depict the kinetic energy 
profiles Et/2t of the corresponding FDS and HDS. 

we have already seen that they markedly influence the 
charge and spin distribution. On the one hand, one sees 
that the AF domains between stripes are characterized 
by a more favorable kinetic energy gain in the SBA, ow- 
ing to a combined effect of a weaker spin polarization and 
larger doubly occupancies than those found in the HFA. 
On the other hand, in the SB approach transverse charge 
fluctuations on the bonds connecting the DWs, with a 
strong reduction of D, and their nearest neighbors, are 
less efficient than the corresponding ones found in the 
HFA, where the fact that DWs are nonmagnetic implies 
much larger D at these sites [cf. Fig. HI a)]. However, in 
the SBA, a small D does not necessarily imply a strong 
reduction of the band narrowing factor. Indeed, the op- 
erator z<ra> describes a sum of the two possible transi- 
tion channels, either of which may accompany an elec- 
tron hopping process^ (singly occupied site) — > (empty 
site), or (doubly occupied site) — > (singly occupied site 
with a time-reversed spin). Therefore, as the hole den- 
sity nh at DW sites is large, the first channel should even 
dominate, which follows from a simple relation, 

e 2 = n h + d 2 . (23) 



that follows for a PM site from the local constraint 10. 
In addition, the HF stripe phase has a twice smaller elec- 
tron density at the DWs than predicted by the SBA. Al- 
together, despite large double occupancies at the DWs, 
the HF kinetic energy gain along the FVS is smaller than 
the SB one, as shown in Fig.JfJb). 

Nevertheless, we have found a confirmation of the 
trend detected in the earlier calculations^*^ performed 
in the HFA — also in the SBA the FVS are mainly stabi- 
lized not by the charge motion along the DWs but rather 
by the hole derealization along the bonds perpendicu- 
lar to the stripes, i.e., by the so-called solitonic mecha- 
nism. This effect is well illustrated by both components, 
Ef and Ef , of the kinetic energy E t , shown in Table 
IIIII Here, one observes a significant anisotropy between 
the HF kinetic energy gain in the x and y-directions, 
SE = Ef — E\. This anisotropy is considerably reduced 
when the local electron correlations are included in the 
SBA, but remains the driving mechanism stabilizing the 
stripes. 

These results should be now compared with those in 
Figs. EJc) andEtd), obtained for the structure with the 
HVS. First of all, note that as in the case of the filled 
stripes, inclusion of the strong electron correlations im- 
proves the kinetic energy gain in the AF domains. How- 
ever, the fact that the half-filled DWs are characterized 
by a smaller hole density than their filled counterparts 
results in an entirely different mechanism being respon- 
sible for their stability. Indeed, Fig. OJd) shows that the 
largest kinetic energy gain is released by the electrons 
propagating along the half-filled DWs, while the trans- 
verse charge fluctuations are less important. This effect 
is particularly strong in the HFA, where one finds a large 
anisotropy between Ef(l x ) and Ef(l x ) at the DWs. We 
ascribe this to the fact that the increase of electron den- 
sity at the nonmagnetic vertical DWs raises double oc- 
cupancy (cf. Figs. |3|andQJ in this approach, and there- 
fore facilitates the on- wall propagation of the electrons. 
Moreover, the nearest neighbor sites of the half-filled DW 
possess a much less quenched S„, as compared to the 
neighbors of the filled one. Therefore, all the electrons 
crossing HVS encounter in this case a stronger on-site 
potential that develops in the AF domains, which acts 
to suppress locally Ef(l x ). A similar amplification (re- 



TABLE III: Comparison of the site- normalized kinetic energy 
contribution for bonds along (10) (Ef) and (01) (Ef) direc- 
tions, as obtained for the phase with either the FVS or HVS 
in the Hubbard model on an 128 x 128 cluster at x = 1/16 
doping and (3t = 1000 within the HFA and SBA. 







FVS 




HVS 


u/t 


method 


Ef/t 


E y t /t 


Ef/t Ef/t 


6 


HFA 


-0.5782 - 


0.5436 


-0.5337 -0.5905 




SBA 


-0.6085 - 


0.5844 


-0.5638 -0.6171 


12 


SBA 


-0.3872 - 


0.3583 


-0.3562 -0.3794 
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FIG. 6: (Color online) Local hole n h (l x ) (top) and magne- 
tization St:(I x ) (second row) density, as well as double oc- 
cupancy D(l x ) (bottom), of the FVS (left) and FDS (right) 
phase found in the Hubbard model on an 128 x 128 cluster for 
doping x = 1/8 at f3t = 1000. Open (filled) circles show the 
HF (SB) results for U = 6t, respectively, while filled squares 
denote the SB data for U = 12t. 



duction) of the on-wall (transverse) kinetic energy gain, 
respectively, we have also found in the SBA. 

We conclude the above analysis by showing in Table llTll 
two kinetic energy components and E\ as found for 
the phase with the HVS. As in the case of the FVS, the 
presence of the HVS introduces a finite anisotropy be- 
tween the kinetic energy gain in the x and y directions 
but while E'f < E\ for the phase with the FVS, the struc- 
ture with the HVS is characterized instead by Ef > Ef . 

Let us now discuss variation of the kinetic energy gain 
shown in Figs, [^fe) and^f) across the unit cell of the 
phase with either the FDS or HDS. Here one finds that 
due to the symmetry E t /2 = Ef — E\ . Therefore, based 
on Fig. EJe) showing substantial kinetic energy gain due 
to transverse hopping processes onto and off the stripe, 
one might expect that such a structure enables the largest 
kinetic energy gain among the stripe phases. In fact, it 
happens only in the strong coupling regime with U = I2t, 



as reported in Table ITT1 In contrast, a weak charge mod- 
ulation induced by the HDS results in a smooth kinetic 
energy profile without a particularly severe gain in the 
vicinity of the DWs. As a consequence, this phase is char- 
acterized by the smallest kinetic energy gain regardless of 
both the approach and the strength of U (see Table ITTjl . 



C. Stripes at intermediate doping x = 1/8 

We turn now to a short discussion of the stripe pro- 
files shown in Figs. |H|and0 obtained at a twice higher 
doping level x = 1/8. The stripe structures at this dop- 
ing confirm the main conclusions of the previous Sec- 
tion concerning the stripe stability. One finds again that 
the transverse (parallel) hopping stabilizes the structures 
with filled (half-filled) vertical DWs, respectively (c/. Fig. 
IBJ. In addition, we point out the most prominent changes 
with respect to the low doping regime. 

First of all, the increased doping acts to shrink the 
distance between the filled DWs according to the for- 
mula d — 1/x, so that the charge (magnetic) unit cell is 
reduced and involves now eight (sixteen) atoms, respec- 
tively, as depicted in Fig. One finds that the charge 
distribution is smoother in the SBA than in the HFA, 
and in the moderate coupling regime U — 6t the highest 
hole density on the nonmagnetic (SV = 0) DW atoms is 
only rih — 0.22 ~ 0.24. This demonstrates that the filled 
stripe phases are truly itinerant systems with modulated 
charge density. 

However, the experimentally established distance be- 
tween stripes in Nd-LSCO at x — 1/8 doping is equal to 
d = 4, 7 implying that de facto less holes are doped per 
one DW site, and consequently DWs are only partly filled 
by holes. Indeed, other local minima are obtained in the 



TABLE IV: 


The same 


as in Table ITU but for x 


= 1/8. 


method 


U/t 


phase 


Eu/t 


Et/t 


F/t 


HFA 


6 


HDS 


0.6321 


-1.2330 


-0.6009 






AF 


0.5148 


-1.1808 


-0.6660 






HVS 


0.5031 


-1.1800 


-0.6769 






FVS 


0.4649 


-1.1777 


-0.7128 






FDS 


0.4432 


-1.1639 


-0.7207 


SBA 


6 


HDS 


0.4815 


-1.2713 


-0.7898 






AF 


0.4563 


-1.2563 


-0.8000 






HVS 


0.4083 


-1.2098 


-0.8014 






FVS 


0.4269 


-1.2470 


-0.8201 






FDS 


0.4123 


-1.2324 


-0.8201 


SBA 


12 


HDS 


0.2130 


-0.7467 


-0.5337 






AF 


0.2726 


-0.8112 


-0.5386 






HVS 


0.2486 


-0.8175 


-0.5689 






FVS 


0.2645 


-0.8396 


-0.5751 






FDS 


0.2588 


-0.8409 


-0.5821 
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FIG. 7: (Color online) The hole density n b (l x ) (top), magne- 
tization density S n (l x ) (second row), and double occupancy 
D(l x ) (bottom) for the half-filled stripe phases: HVS (left) 
and HDS (right). Open (filled) circles show the HF (SB) re- 
sults for U — 6t, respectively, while filled squares denote the 
SB data for U = 12t. Doping, cluster size and temperature 
as in Fig. HJ 

SB calculations for the charge (magnetic) unit cells of the 
half-filled stripe phases, containing four (eight) atoms, 
respectively, as shown in Fig. Also these phases are 
itinerant, with the SB hole densities at U = 6t less than 
0.2 at the DW sites for both half-filled stripe phases. 
Moreover, correlations implemented in the SB approach 
guarantee that electrons avoid each other and the double 
occupancy is the lowest D ~ 0.07 at the half-filled verti- 
cal DW sites, while in the HFA one finds instead a large 
value Dq ~ 0.13. However, it is remarkable that the SB 
double occupancy D is practically site independent in the 
phase with the HVS, in spite of rather different electron 
densities at the DW sites and within the AF domains. 

Next, we address a generic crossover found in the SBA 
from the FDS to their vertical counterparts taking place 
upon doping. We have seen from the respective charge 
distributions in Fig.|g|that a certain amount of holes pro- 
liferates out of the DWs, reducing double occupancies in 
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FIG. 8: (Color online) Kinetic energy Ef(l x ) and E%(l x ) pro- 
jected on the bonds along the x and y direction, as obtained 
for the phase with either the FVS (a,b), or HVS (c,d), shown 
in Figs. and respectively. Open (filled) circles show the 
HF (SB) results for U = 6t, respectively, while filled squares 
denotes the SB data for U = 12i. 



the AF domains. Consequently, the on-site energy gains 
become less important, which explains the transition to- 
wards the phase with the FVS, being a more favorable 
structure for the charge dynamics. In fact, for U = 6t, 
this crossover occurs precisely at x = 1/8, so that both 
phases become degenerate, as reported in Table llVl while 
in the strong coupling regime U = 12t one recovers in 
the ground state the phase with the FDS. Moreover, as 
the doped holes are now only weakly localized at the 
DW sites, a nonuniform charge distribution in the AF 
domains amplifies the anisotropy between the kinetic en- 
ergy gain Ef and Ef (cf. Table©. 



D. Electronic structure of stripe phases 

Let us analyze in more detail the electronic structure 
obtained for the stripe phases in the low doping regime. 
Since the stripe phases consist of a superposition of mag- 
netic domains with a 2D character, and domain walls of 
ID character, it is worth determining to which part of the 



TABLE V: The same as in Table |m| but for x = 1/8. 







FVS 




HVS 


u/t 


method 


Ef/t 


EUt 


Ef/t E y Jt 


6 


HFA 


-0.6234 - 


0.5543 


-0.5331 -0.6469 




SBA 


-0.6468 - 


0.6002 


-0.5504 -0.6593 


12 


SBA 


-0.4486 - 


0.3910 


-0.3848 -0.4327 
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FIG. 9: Density of states N(u)) of the phase with: (a,b) HVS, 
shown in the left panels of Fig.Q] (c,d) FVS, depicted in the 
left panels of Fig. El as obtained within the HFA (left) and 
the SBA (right) in the moderate coupling regime U — 6t for 
doping x = 1/16. 
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FIG. 10: Density of states N(uj) of the phase with: (a,b) 
HVS, shown in the left panels of Fig. (c,d) FVS, depicted 
in the left panels of Fig. HO as obtained within the HFA (left) 
and SBA (right) in the moderate coupling regime U — 6t for 
doping x = 1/8. 



spectrum each of them contribute, and to which extend 
they mix. The DOS of a pure 2D AF phase exhibits a 
gap (A ~ U) separating two van Hove singularities lying 
at the top (bottom) of the lower (upper) Hubbard sub- 
bands, here representing the dynamical Hubbard bands. 
In ID case the van Hove singularities are located at the 
lower and upper band edges. 

Let us consider the actual DOS for the phase with the 
HVS in the moderate coupling regime U — 6t — it con- 
sists indeed of two Hubbard subbands at low and high 
energies, lower Hubbard band (LHB) and upper Hub- 
bard band (UHB), separated by a gap of ~ 5t in the 
HFA and ~ 2.6f in the SBA, see Figs.EJa) andEfb). So 
strong reduction of the gap with a similar value of the 
magnetization in the AF domains is the first correlation 
effect identified in the electronic structure. 

Furthermore, the formation of stripes with nonmag- 
netic sites results in additional electronic states within 
the Mott-Hubbard gap. These new structures contain 
parts of the spectral density below the Fermi energy, both 
in the HFA and in the SBA, indicating that half-filled 
stripe phases are metallic. Note also that the nonmag- 
netic stripes hybridize only weakly with the neighboring 
AF domains, since they are mainly stabilized by the elec- 
trons propagating along them, so this part of N(oj) is 
reminiscent of the tight binding DOS for a linear chain 
with nearest neighbor hopping, and has the character- 
istic peaks at the edges. We recall that in ID case the 
van Hove singularities are located at the lower and upper 
band edges, as found here for the in-gap bands. Similar 
to the Hubbard subbands, also the partly filled electronic 
subbands in the gap are substantially narrowed by the 



correlation effects in the SBA. 

In the case of filled vertical stripe phases shown in Figs. 
OJc) and OJd) , the states in the gap are empty. This 
was actually pointed out rather early as the microscopic 
reason responsible for the stability of such phases in the 
HFA4i The most prominent feature of the DOS for FVS 
phase is again a striking similarity of the mid-gap part 
of N(lu) to the DOS of a linear chain. The width is 
somewhat reduced indicating a more localized character 
of these states. In contrast with the phase with the HVS, 
a stronger hybridization with the AF domains influences 
the UHB so that a certain amount of its states is shifted 
towards a lower energy and forms a separate segment of 
N(u) [cf. Fig. He)]. " 

Finally, in Fig. Eld) we show the DOS of the phase 
with the FVS which follows from the SBA. In this case 
the charge and spin density profiles are more spread out 
and a weaker anisotropy between the tranverse and paral- 
lel electron motion, as compared to the HF ones, results 
in a less clear character of the mid-gap part of N(u>), 
with a large peak in the middle being now a characteris- 
tic feature of the 2D tight binding DOS. Notice also an 
insulating nature of both filled stripes. Indeed, as the 
mid-gap states are now entirely empty, the Fermi energy 
lies inside the gap between the highest occupied state of 
the LHB and the bottom of the mid-gap band. 

For the higher doping x = 1/8 one finds similar fea- 
tures in the DOS of both HVS and FVS phases. The 
spectral intensity is here redistributed, showing a sys- 
tematic transfer of the spectral weight from the UHB 
into these states and the LHB upon doping, in agree- 
ment with a strong coupling perturbation theory for the 
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Hubbard model^ In addition, the spectral weight is also 
transfered into the mid-gap part of N(ui), which grows 
with an increasing number of stripes in the cluster. 



IV. STRIPE PHASES IN THE EXTENDED 
HOPPING MODEL 

A. Crossover to vertical stripe phases 

Thus far we have worked out that in the strongly cor- 
related regime relevant for the cuprates a structure with 
the FDS (HVS) is the lowest in energy among the stripe 
structures with the filled (half-filled) DWs, respectively. 
Therefore, an interesting question occurs — which mi- 
croscopic parameters decide whether the phase with the 
FDS or the one with the HVS is more stable. To clarify 
this point we have investigated in the SBA the competi- 
tion between the stripe phases in the t-t'-U model with 
U = 12t. The chosen value of the Coulomb interaction 
gives the ratio J/t = 1/3 which corresponds to the phys- 
ical value in the cuprates 

Fig- EI a) shows the effect of increasing next-nearest 
neighbor hopping \t'\ on the free energy F of all four 
stripe phases, as well as of the uniform AF phase at the 
hole doping x = 1/16 and U = 12t. Here, the most 
striking result is that while increasing \t'\ only weakly 
influences the energy of the AF phase, it modifies sub- 
stantially F of the stripe phases. Moreover, it clearly 
stabilizes half-filled stripes reducing simultaneously sta- 
bility of the filled DWs. To appreciate better the en- 
hanced stability of the half-filled stripes with respect to 
the latter ones, we analyze the average kinetic energy 
Ef per diagonal bond as a function of increasing \t'\, ob- 
tained in the AF and PM phases for a few doping levels: 
x = 1/16 (dot-dashed line), x = 1/8 (dashed line), and 
x = 1/5 (solid line) (see Fig. Il2|l . One finds that a nega- 
tive next-nearest neighbor hopping (t'/t < 0) yields first 
a positive kinetic energy contribution with the loss being 
proportional to doping. On the one hand, hopping pro- 
cesses associated with t' are partly suppressed in the AF 
phase by both the reduction of double occupancies and 
the increase of spin polarization. On the other hand, in 
the PM phase, these processes might be optimized only 
by the first effect, which results in a much more severe 
loss of the energy. However, these trends also involve 
a reduction of the nearest neighbor hopping processes 
and therefore they cannot be maintained above a certain 
critical value of t' . As a result, Ef changes its sign and 
becomes negative (cf. Fig. I12JI . 

Turning back to the analysis of the stripe stability, the 
fact that the energy loss due to the t' term follows pre- 
dominately from the hopping processes at the nonmag- 
netic DWs explains immediately why a finite t' promotes 
partially filled stripes. Indeed, as the average hole den- 
sity at the half-filled DWs is smaller than at the filled 
ones (cf. Figs. El and SJ, one should lose less kinetic en- 
ergy in the former case, especially in the phase with the 
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FIG. 11: (Color online) Site-normalized SB free energy F 
of the phase with the FVS, FDS, HVS, and HDS, as well 
as of the AF phase, as a function of increasing next-nearest 
neighbor hopping [t | for U = 12t in the: (a) underdoped 
regime x = 1/16; (b) moderately doped regime x = 1/8. 



HDS, where the reduction of the hole density at the DWs 
is the strongest. Consequently, as depicted in Fig. HTT a,). 
this phase is characterized by the largest free energy gain 
upon increasing \t'\. Further, increasing \t'\ strongly in- 
terferes with the solitonic mechanism-A stabilizing the 
FVS by the transverse hopping and hence such a struc- 
ture becomes unstable at t'/t ~ —0.248. In contrast, nar- 
rower FDS are stabilized to a lesser degree by this mech- 
anism remaining a stable solution in the entire region 
< \t'/t\ < 0.4. Finally, the energy difference between 
the phase with the HVS and the one with the FDS grad- 
ually diminishes so that the former becomes the ground 
state of the system at t'/t ~ —0.16. 

In order to get a more quantitative insight into the 
mechanism of the crossover between the above phases, we 
show in Table Ivfl decomposition of their free energy into 
the on-site Coulomb energy En as well as into both ki- 
netic energy contributions E t and E t > . On the one hand, 
the robust stability of the FDS in the absence oft', follows 
mainly from a large kinetic energy gain E t (cf. Table ITTll. 
In contrast, at the expense of E t , the structure with the 
HVS better optimizes the double occupancy energy Ejj, 
being however less stable than the one with the FDS. 
On the other hand, as we have already found, a negative 
next-nearest neighbor hopping tends to yield a positive 
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FIG. 12: (Color online) Average kinetic energy E t i per diag- 
onal bond as a function of increasing next-nearest neighbor 
hopping \t'\ for: (a) AF phase, and (b) PM phase, as obtained 
for U — 12t. Dot-dashed, dashed, and solid line shows the re- 
sults for x = 1/16, x — 1/8, and x = 1/5 doping, respectively. 



FIG. 13: Density of states JV(lj) of the ground state with: 
(a) FDS for x = 1/16 and t' = 0; (b) HVS for x = 1/8 
and t' = — 0.3t, as obtained for large Coulomb interaction 
U = I2t. 



kinetic energy contribution, but the actual sign depends 
both on the value of t' and on the filling of the DWs. In- 
deed, the presence of the half-filled stripes clearly helps to 
alterate the disadvantageous sign of E t > so that it is nega- 
tive already for t' = — 0.3i (cf. Table E} . Conversely, in 
the phase with the FDS, in spite of a stronger reduction 
of double occupancies with a concomitant detriment of 
Et, derealization of electrons due to a finite t' = —0.3t 
still costs a small amount of energy. All these features 
act to stabilize the HVS in the ground state. 

Based on the above discussion one could expect that 
at a higher doping, the HVS would take over the FDS 
for a much smaller value of \t'\. Indeed, as shown in 
Table EH x = 1/8 provides a more significant energy 
lowering of the phase with the HVS due to a twice larger 
Ef gain. Simultaneously, the energy loss in the structure 
with the FDS due to the t' processes is twice as large as 
for x — 1/16. Consequently, increasing \t'\ destabilizes 
now easier the FVS and one finally loses this solution 
already at t'/t ~ —0.247, as illustrated in Fig. fTTf bl. 
Nevertheless, one reads off from Fig. lllf b) that a change 
of the ground state upon increasing \t'\ takes place in 
this case at even slightly larger value \t'\ ~ 0.17i. The 
explanation is contained in Table IVII a higher doping 
diminishes the double occupancies and additionally it is 
reflected in a higher mobility of the electrons. As a con- 



TABLE VI: Free energy F per site, on-site Coulomb en- 
ergy Eu as well as kinetic energy Et and E[ contributions, as 
obtained for the phase with either the HVS or FDS in the ex- 
tended Hubbard model from the SBA on an 128 x 128 cluster. 
Parameters: U = 12t, t' = -Q2t and fit = 1000. 



X 


phase 


Eu/t 


Et/t 


E t ,/t 


Fit 


1/16 


FDS 


0.2775 


-0.7292 


0.0039 


-0.4478 




HVS 


0.2767 


-0.7225 


-0.0071 


-0.4529 


1/8 


FDS 


0.2332 


-0.8064 


0.0077 


-0.5655 




HVS 


0.2300 


-0.7911 


-0.0138 


-0.5749 



sequence, the relative energy 

SF = Fpvs — -FkvS) (24) 

varies, in the absence of t' , from SF = — 0.0065£ for 
x = 1/16 to SF = -0.0132i at x = 1/8, which demon- 
strates an increasing stability of the FDS. Therefore, the 
condition of having the crossover at x = 1/8 requires 
overcompensation of the larger |<5F| and thus the transi- 
tion is now slightly moved towards the larger \t'\. 

Altogether, Fig. El enables scenario in which doping- 
induced increase in the amplitude of t' results in a drastic 
change in the spin modulation from the diagonal to verti- 
cal/horizontal one, as observed experimentally in LSCO 
at x ~ 0.06. The conjecture that \t'\ is suppressed in 
underdoped regime and grows with increasing x, seems 
to be also supported by the recent ED studies showing 
that t' helps to stabilize the superconductivity in the op- 
timally and overdoped regime but it is harmful in the 
underdoped regime4ft Moreover, inelastic neutron scat- 
tering experiments on YBCO, a bilayered compound with 
a larger than LSCO next-nearest neighbor hopping term 
\t'\ ~ 0.3i(2& have established the presence of the IC verti- 
cal/horizontal spin fluctuations throughout its entire su- 
perconducting regimejAA implying that such a modulation 
is: (i) realized easier for larger \t'\, and (ii) more advan- 
tageous for the superconductivity than the diagonal one. 
Note, however, that it remain puzzling why the average 
filling of diagonal DWs changes from one to (l/y/2) hole 
per one DW site in the very low doping regime x ~ 0.03. 
Unfortunately, we could not address this issue here as it 
requires calculations for still larger unit cells which could 
not be performed at present. 



B. Changes in the electronic structure due to 
finite t 

Remarkably, in the present approach, the two possible 
ground states have entirely different physical properties 
being either insulating or metallic. Indeed, the DOS of 
the ground state with the FDS obtained in the Hubbard 
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FIG. 14: (Color online) Electronic structure of the phase with 
the HVS as a function of parallel momentum k y , calculated 
from the unit cell shown in Fig. ffl^a) on an 8x8 cluster within 
the SB method for U/t = 12, and: (a) t' = 0, (b) t'/t = -0.3. 
Solid (dashed) lines correspond to the bulk (mid-gap) bands, 
respectively, whereas the dotted lines indicate the Fermi level. 



model with t' = at x = 1/16 consists of clearly seen two 
distinct maxima corresponding to the Hubbard bands, 
while the Fermi energy falls within a gap which results 
from the magnetic order [cf. Fig. IT^T a)]. Interestingly, 
due to a strong Coulomb repulsion U = Y2t and con- 
sequently small band narrowing factors, instead of the 
ones as for U — 6t, one finds here four separated midgap 
bands. In contrast, as shown in Fig. IT^ h). the Fermi 
energy crosses the midgap states of the structure with 
the HVS, which is the ground state of the Hamiltonian 
with x = 1/8 and t' = — 0.3t, enabling charge transport 
in agreement with the data for LSCO. 

We also note that, in agreement with Ref. l4ll a nega- 
tive t' leads to a distinct broadening of the partially filled 
mid-gap band induced by the HVS and consequently 
it shifts these states to a lower energy, as depicted in 
Fig. Conversely, the energy gain due to such a mod- 
ification of the electronic structure is not possible in the 
case of the FDS as their mid-gap states are entirely un- 
occupied. Thus, the puzzling role of t' in promoting the 
half-filled DWs is clarified. 



V. DISCUSSION AND SUMMARY 

We have developed a simple but powerful SB approach 
which allows one to investigate various stripe phases with 
a large unit cell and carry out the calculations on a large 
(128 x 128) cluster. Having presented the theoretical 
framework in Sec. [HJ we have shown in Sec. Ill Dl that 
our method provides a unique opportunity to obtain un- 
biased results at low temperature (it = 1000, as well as to 
eliminate the role of finite size effects. The electron cor- 
relation effects beyond the HF were implemented within 
the SB method. Therefore, the stripe phases found in 
the present approach are stabilized not due to particular 
boundary conditions, but they represent a generic insta- 



bility of the strongly correlated electron system doped by 
holes. 

Then, in Sec. lIIII we have analyzed the stability of the 
idealized filled as well as the half-filled stripes in the Hub- 
bard model at two representative doping levels x = 1/16 
and x = 1/8. However, the true ground state could corre- 
spond to neither filled nor half-filled stripes, as the opti- 
mal filling might vary with doping. By comparing the SB 
charge and spin density profiles with the ones obtained 
in the HFA in moderate coupling regime with U = 6t, we 
have emphasized the role of a proper treatment of elec- 
tron correlations. In particular, we have found that the 
largest correction of the HF value is obtained at the DW 
unpolarized sites, where the double occupancy D shows 
distinct minima, which allows one to optimize the on-site 
Coulomb energy Ejj even without a great reduction of 
the actual electron density. Simultaneously, an enhanced 
value of D in SB approach suppresses partially the spin 
polarization of the atoms within the AF domains, which 
enables more intersite excitations and leads to a more 
favorable kinetic energy gain. Taken together, these two 
effects are responsible for smoother SB charge and spin 
density profiles with respect to the ones found in the 
HFA. Moreover, we have shown that at strong coupling 
with U = Ylt the SB double occupancy profile is interpo- 
lating between 4t 2 /U 2 deep in the magnetic domains and 
the value expected in the PM phase at the hole density 
obtained at the DW sites. 

However, the most prominent result of our extensive 
studies of the kinetic energy was a demonstration of a 
close relationship between the filling and the direction of 
the largest kinetic energy gain. Indeed, we have shown 
that while deep within the magnetic domains the kinetic 
energy in both directions is given by — At 2 /U, as expected 
from large U expansions^ in the neighborhood of DWs 
the physical situation is much more complex and depends 
on the microscopic properties of the considered stripe 
phase. Especially, it turns out that the HVS are pre- 
dominately stabilized by the electrons propagating along 
them, whereas their filled counterparts might be consid- 
ered as a set of weakly coupled solitons as their stabil- 
ity follows to a large extent from the transverse hopping 
across them. 

Finally, proceeding from the experimental motivation, 
based on the doping-induced change of the spatial orien- 
tation of the DWs in the cuprates, we have investigated 
the effect of the effective next-nearest neighbor hopping t' 
on the stripe stability. Remarkably, we have found that a 
subtle increase of the amplitude of t' with increasing dop- 
ing could explain the observed transition from the filled 
insulating diagonal stripes towards the half-filled metallic 
vertical/horizontal ones in the cuprates. As the effective 
hopping t' is different for hole-doped and electron-doped 
cuprates and the doped holes delocalize the magnetic 
moments at Cu sites, it might be expected that \t'\ in- 
creases with doping. This poses an interesting physi- 
cal problem for the parameters of the effective one-band 
model for the cuprates which should be resolved by future 
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studies. Such a doping dependent t' may also result from 
the change of lattice parameters upon doping observed 
in most superconducting cupratesi5£ 
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